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|- We propose a theory of the dynamics of polymers in dilute solution, in which the popular Zimm 

D ' and Rouse models are just limiting cases of infinitely large and small draining parameter. The 

equation of motion for the polymer segments (beads) is solved together with Brinkman's equation 
for the solvent velocity that takes into account the presence of other polymer coils in the solution. 
The equation for the polymer normal modes is obtained and the relevant time correlation functions 
are found. A tendency to the time-dependent hydrodynamic screening is demonstrated on the 
diffusion of the polymers as well as on the relaxation of their internal modes. With the growing 
concentration of the coils in solution they both show a transition to the exactly Rouse behavior. 
The shear viscosity of the solution, the Huggins coefficient and other quantities are calculated and 
CZ) , shown to be notably different from the known results. 



in 



o 
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Interest in dilute polymer solutions arises primarily from their importance in the characterization of polymers, their 
interaction with solvent and from a fundamental interest in understanding macromolecular response to hydrodynamic 
forces, free from the complications of intermolecular entanglements. Despite five decades of investigations, these 
questions are not entirely understood. So, there is a systematic discrepancy between the dynamic scattering data and 
the theory [HE El- The existing theories give different results for the viscosity of dilute polymer solutions (see, e.g., 
. 0> § and the citations there), the observed monomer motion in single polymer chains cannot be explained by the 
available theories , the time dependence of the hydrodynamic screening in solutions is not explained 7] , etc. For 
other problems we refer the reader to the recent review The aim of this work was to contribute to the solution of 
. some of these problems by developing a phenomenological bead-spring theory of the diffusion of an individual "test" 
polymer and the relaxation of its internal modes in solution of unentangled polymers. Our approach differs from 
the traditional ones in several main points. First, the joint Rouse-Zimm theory is exploited |9j. In the literature, 
following de Gennes 01 j it is assumed that at 9 conditions the Zimm modes with the dispersion of relaxation times 
T pZ ~ p~ 3 ^ 2 , where p is the mode number, at low frequencies always must dominate the Rouse modes with t p r ~ p~ 2 ; 
we show that both these contributions to the polymer characteristics must be taken into account. Next, the internal 
modes are distributed discretely (the assumption of their continuous distribution with respect to p is true only in a 
restricted time domain and often leads to incorrect interpretation of experimental data The formalism from 

our theory with the time-dependent hydrodynamics of polymers has been adopted 0, IT^ | and, finally, the Brinkman's 
theory [lj] for the flow in porous media is used to take into account the influence of other coils on the studied 
test polymer. The presented theory has the following limitations. The considered time scales are t >> tr = R 2 p/r], 
, where R is the hydrodynamic radius of the polymer, p is the density and r\ the viscosity of the solvent. This means 
5_i ■ that the effects of hydrodynamic memory (or the viscous aftereffect) are neg lected HUJ,UJ. The distribution of the 
coils in solution is considered to be stationary (this is justified at least for the times t << T£>; the choice of this time 
scale is possible since always r p << Try, where to = R 2 /D is the characteristic time of the coil diffusion with the 
diffusion coefficient D). Our theory is also restricted to solvents 0; generalizations to other cases require knowledge 
of the equilibrium distribution of the segments when the exclude volume interactions are taken into account. As 
already mentioned, only solutions of unentangled polymers are considered. We are thus limited to concentrations of 
the chains c < 1 / [77] , where [rj\ is the intrinsic viscosity and c the number of polymers per unit volume 8J. In spite 
of all these restrictions and other ones like those that we do not consider the internal viscosity of polymers and the 
self-entanglements (the importance and even the reality of these interactions are uncertain |8j), we believe that the 
results of our theory could be of interest. In particular, we have found new expressions for the quantities describing 
the behavior of flexible polymers in solution, such as the diffusion coefficient of the coil, the relaxation times of the 
internal modes, the viscosity of the solution, and the Huggins coefficient. These quantities have been obtained from a 
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generalized Rouse-Zimm equation for the position vectors of the polymer segments and the Oseen tensor describing the 
velocity field of the solvent due to perturbation. Finally, our theory describes in a simple manner the hydrodynamic 
screening, i.e. the concentration and time-dependent transition between the Zimm and (as distinct from the previous 
theories) exact Rouse behavior of the polymer. 



II. THE ROUSE-ZIMM-BRINKMAN THEORY OF POLYMER DYNAMICS 

We choose one coil as a "test" polymer. The equation of motion of its nth segment is 

M%^7t + 7f + 7V (i) 

Here, ~a?„ is the position vector of the segment (a spherical bead) from the N ones constituting the polymer, M is 
the bead mass, ~f n is the force from the neighboring beads along the chain, / n the random force due to the motion 
of the molecules of solvent, and ~f n is the Stokes friction force on the bead during its motion in the solvent P.ll6|: 



n v , V, v 

v[xn) 



(2) 



where if denotes the velocity of the solvent in the place of the nth bead due to the motion of other beads. The friction 
coefficient on the bead with radius b is £ = 67rn6. This expression holds in the case of steady flow and takes into 
account the hydrodynamic interaction. In a more general case with the hydrodynamic memory 0, IT2I Il5| the force 
(2) should be replaced by the Boussinesq force and equation (1) has to be solved together with the nonstationary 
hydrodynamic equations for the macroscopic velocity of the solvent. To take into account the presence of other 
polymers in solution, we use the Brinkman's work |l3j | (see also in which a polymer is considered as a porous 
medium. In our approach all the solution is such a medium with coils being obstacles to the solvent flow. Then in the 
right hand side of the Navier-Stokes equation a term —K 2 rf%f has to be added, where k~ 2 is the solvent permeability. 
This term has a sense of the average value of the force acting on the liquid in an element of volume dV, provided the 
average number of polymers in solution per unit volume is c; then K 2 n = cf, where / is the friction factor on one coil. 
Thus, for an incompressible solvent (S/lf = 0)we have to solve the equation 

p—^r— = — Vp + rjAlf — n rfvf + ~0 . (3) 
at 

Here p is the pressure and Tfi. is the density of the force from the beads of the studied polymer on the solvent , 

= -£7i[ r (4) 

n 

To solve this equation is a difficult problem since the polymer chains are mobile. However, restricting ourselves to the 
times much shorter than td, the concentration c can be assumed constant. The above equations then describe the 
motion of one bead in the solvent with the effective influence of other coils on the motion of the solvent flow. This 
problem can be transformed to that solved already in 01 ( see a l so El El)- The velocity field can be in the Fourier 
representation in the time written as follows: 

< (f*) = / dT" J2 HSp - fp • (5) 

Here the analog of the Oseen tensor is 

K/3 (^) - ASafi + Br a r p r- 2 , (6) 
A = {Smjr)- 1 \e- y - y [(l - e~ v ) 2/" 1 ]"} , 

B = (STrnr)- 1 {e"" + 3y [(l - e~») 2T 1 ]"} , (7) 

V = r Xi X 2 = k2 ~ i^Q/Vi an d the prime means the differentiation with respect to y. In theparticular case to = 
and for permeable solvent, k = 0, equation (6) coincides with the well-known result of Zimm Using this solution, 
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a generalization of the Rouse-Zimm equation can be obtained from the equation of motion 
the Oseen tensor over the equilibrium Gaussian distribution of the beads 0, ^| gives 



( H af3nm) = S a ph U (fl - 77l) , ~t n m = X 



The preaveraging of 



(8) 



{n-m) = (67r 3 |n-m|) 1/2 (rja^ 1 [l - spnz exp (z 2 ) erfc(z)] . 

1 /2 

Here a is the mean square distance between the beads along the chain and z = x a (l n ~~ m 

|/6) 7 . Then in the 

continuum approximation with respect to the variable n the new Rouse-Zimm equation reads 



1 



3k B T d 2 ^ \n) 



<9n 2 



(9) 



JV 



+ / dmh^fa — m) 



3k B Td 2 lt u (m) 



dm 2 



+ Mw 2 ^(r7i) + 7"(m) 



It is solved with the help of the Fourier transformation (FT) in n, taking into account the boundary conditions at the 
ends of the chain n = at n = 0, N: !& w {n) — ~y*$ + 2 X)p>i ~l?p cos(irnp/N). The inverse FT then yields 



where 



5£ = £ [1 + (2 - M N£h% p ] \ K p = 3k B T 
and the Oseen matrix is 0, 

1 



1 + Xp 



iV 



Ttriay/MNp l + (l + Xp y 



O X°' P = l' 2 ' • 

37rp 



(10) 



(11) 



(12) 
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fl - exp x^erfcx^ 



Xa. 



Using the fluctuation-dissipation theorem, the time correlation functions of the normal modes are 

k B T r°° . ReS^ 



il> P (t) = (y ap (0)y ap (t)) 



(2 - Spo) ttN J_ 



dui COS Ujt- 



-iuiSf - Muj 2 + K P \ 



(13) 



(14) 



A. Diffusion of the coil 



In the stationary limit u> = so that x — K - Then the preaveraged Oseen tensor (6) is 

/riw \ _ / exp(-xr) \ 



(15) 



The quantity 1/k can be thus (for small nr only) considered as a screening length. For an individual polymer we had 
(p = in equation (14)) 



M0)-Mt) = Dt 



(16) 



with the diffusion coefficient D = Dr + Dz (R and Z stay for the well-known Rouse and Zimm limits Q. Now instead 
of equation (12) we have /iq with xo = kRq (Rq being the gyration radius), the diffusion coefficient depends on the 
concentration of the coils c, 



D = D R + D z (c), 



(17) 
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(Dz(0) = Dz) and consists of the Rouse (independent on the presence of other polymers) and the Zimm contributions. 
The latter one can be expressed in the form 



Dz(c) 

where /(c) is a "universal" function for every polymer: 



Dzf(c) 



4xo 



1 - 



V^Xo xl 



(l - expxoerfcxo 



(18) 



(19) 



The dependence of the permeability on the concentration is estimated as follows. The friction coefficient in the 
quantity k 2 — cf /rj from equation (3) can be determined using the Einstein relation D = ksT/ f. In such a picture 



K 2 = 



27v^r c 
16 R 2 G 



4V2h 



(20) 



Then the values of k and xo depend on the draining parameter h — 2(3N/Tr) 1 / 2 b/a (if h » 1, the dynamics is of 
the Zimm type, for h << 1 we deal with the Rouse polymers). The quantity c = denotes the number of 

polymers per the volume of a sphere with the radius Rg- With the increase of c the Zimm term decreases and for 
large c (small permeability k when \o >> 1) it becomes ~ I/v^j 



Dz(c) 



2k B T 1 



nr/Na 2 n 

The realistic case of small c corresponds to xo = kRq << 1 when 

D z (c) = k B Th° Q0 {c) - D z (l 



8^ 



kR g 



(21) 



(22) 



The concentration dependent correction to Dz is thus proportional to \fc and differs from other results (compare 
Ref. |l8j and citations there, where this correction is ~ c). The behavior of a free polymer depends on the draining 
parameter h. If h is large, the Zimm polymer (at c — 0) with growing c should change its behavior to the diffusion 
with the Rouse coefficient Dr. 



B. Dynamics of internal modes 

In the stationary case (oj = 0) and at zero concentration (k — 0) the diagonal elements of the Oseen matrix are well 
known 4]. Now h^ p from equation (13) depend on c. The internal modes relax exponentially as in previous theories, 
ip p (t) oc exp(— \t\/T p ), but their relaxation rates consist of the Rouse contribution and the concentration-dependent 
Zimm part, 

1 1 1 

(23) 



Tp(c) T pR Tpz(c) 

where t p r and t p z(0) = r p z are given in Refs. [H Il6j and 



X P 

which behaves as 

N 



T M i i + (i + x P ) 2 _ 

TPZ(C) = 2 l + y„ TpZ > (24) 



T P z{c) = T pZ ( 1 + ^kV-... ) (25) 



if c — ► 0, and (although unrealistic), as c — > oo one has 



I (No 2 ) 2 X] 

r ^Y^ = W (26) 



Note that for the internal modes the draining parameter depends on the mode number p: h(p) = t p r/t p z — h/y/p. 
The "universal" dependence of t p z{c)/t p z{Q) on \v (24) indicates that with the growing c every polymer shows a 
tendency to become the Rouse one. 
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C. Steady state viscosity and the Huggins coefficient 

The shear viscosity of the solution can be calculated from the formula 0, 01 Il6| 

V(c) =V + -jk B Tc^T p {c). 

P =i 



(27) 



Using equation (24), in the Rouse limit we have the familiar result |j| 77(c) — 77 = irN 2 a 2 bcr]/6. In the Zimm limit at 
small concentrations 



= '(VNa)C - 
2^ \ ' \2 



(28) 



0.425c (iVa 2 ) 3/2 1 + 0. 140c (Na 



,3/2 



where i?z is the Zimm hydrodynamic radius jj] and C is the Riemann zeta function. The first term coincides with 
the known result Q. A more general expression for the viscosity, following from equations (27) and (24), is 



77(c) 



— = -N 2 a 2 bc\^ I 1 



1 



2/i 



1 



Xp 



vpi + (i+x P r 



At very low concentrations when x p << 1 one has 



(29) 



(30) 



Due to the dependence on /i the difference between this and the classical result can be notable. So, for a polymer 
with small h the ratio of the intrinsic viscosity [77]/! = lim c ^o[ 7 ?( c )~ ? 7]/( 7 ?c) at h < 1 (when the polymer is assumed to be 
the Rouse one) to that with h = changes as a function of h from 1 to w 0.55, at h = 0.5 being 30 per cent smaller than 
in the case of a pure Rouse polymer. For a very large h the intrinsic viscosity is [?/]h>>i = 3y^2/irRQ^(3/2) — 6.253Rq. 
Considering the viscosity normalized to this expression, one can find that even for rather large h the difference from 
the traditional result for the pure Zimm polymer is significant. So, at h = 10 it represents some 25% and it is still 
about 10% even for h such large as 50. 

One of the important rheological parameters of polymer solutions is the Huggins coefficient kn ■ It can be determined 
from the general expression for the viscosity (29), using the intrinsic viscosity [77] ^ at zero concentration (see equation 
(30)): 



We find 



3tt 

23/2 



oc ^ 

E ^772 



p=l 



rjc 



h 



[rj\ h (1 + k H [rf\hC 



AV2h 



For large h (the Zimm case) one thus has 



k 



H 



3tt 

23/2 



C" 



00 1 



0.655. 



h 



(31) 



(32) 



(33) 



This value differs from the literature results, see, e.g., 3 where ku = 0.757 is given; in Ref. fr^l one finds kn = 0.6949, 
and in |2(j the calculations gave the value 0.3787. The Freed and Edwards theory pH l22fl possesses an intrinsic 
viscosity, which is inconsistent with the Kirkwood-Riseman steady-state limit and gives the hydrodynamic screening 
even for infinitely dilute solutions (the discussion of this question has been given already in the paper [f|). 

As h — > (the Rouse case), kn approaches zero as kn ~ 27rft,£(7/2)£~2(2) and when h grows, the Huggins coefficient 
slowly converges to the Zimm limit (33). The difference from this limit is significant in a broad region of h, e.g., with 



the maximum ss 1.27 of the function kn jkuzimm at /i = 3, and with kulkuzi 



1.15 for h = 20. 
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D. Monomer motion 



In connection with the unresolved problem of the dynamic nature of hydrodynamic screening in polymer solutions 
(see Introduction), it is of special interest to consider the time-dependent quantities describing the polymer behavior. 
Among such quantities, the relaxation modulus, which determines the shear stress at shear flows can be easily studied 
since it is given simply by a sum of exponentials containing the relaxation times from equation (23) [J, |8( . Here we 
shall briefly focus on simplest (but observable |||) motion of the end monomer within a polymer coil and calculate its 
mean square displacement (MSD). The MSD part due to internal modes is P. ITll] 



, 2 . AN a 2 ^ 1 



9 

^1 P 



1 — exp 



Tp(c) 



(34) 



As already shown, with growing concentration c every polymer tends to behave as a Rouse one, which is due to 
the decrease of the Zimm contribution to the relaxation rates t~ . The time dependence of this screening is well 
displayed considering, e.g., the ratio of the Rouse part of the MSD (i.e. that if the polymer was the pure Rouse 
one, h — 0) to the total MSD in the joint Rouse-Zimm model. This function, (r 2 (t))i n t .r/ {r 2 (t))i n t, depends on 
the draining parameter h, the concentration c, and the time. With the growing t the above relation converges to 
unity showing the transition to the Rouse behavior. For example, at a concentration c = 0.1 and h = 10 we have 
(r 2 (t))i nt} R / (r 2 (t))int ~ 0.75 at t ~ t\r, at t — 2tir the difference from the Rouse MSD is only about 10%, and 
at t = 5tir the initially Zimm polymer becomes indistinguishable from the Rouse one. When the same relation is 
considered as a function of c for different times, one sees that the tendency to approach the Rouse limit with the 
increase of c is more and more expressed as the time growths. At long times, as expected, the polymer behaves as 
the Rouse one already at small concentrations. 



III. CONCLUSION 



The behavior of complex polymer systems that are attractive due to their unusual properties and numerous appli- 
cations cannot be understood without understanding the behavior of a single polymer in a liquid, its interaction with 
the solvent and with other polymers in dilute solutions. Even in situations when we deal with solutions of flexible and 
unentangled polymers, a number of open questions exists for years and new "puzzles" appear. To our opinion, some of 
the problems are only due to inappropriate use of the existing theories and to a great influence of simple " universal" 
laws of polymer behavior, such as the famous fc 3 law for the first cumulant of the dynamic structure factor or the t a 
laws for the monomer MSD, where a = 1/2 for the Rouse polymer and 2/3 for the Zimm one. A closer look at these 
laws shows, however, that their application to real situations is rather restricted and often they do not correspond to 
experimental conditions. The model developed in the present work is not particularly new. Partially (when dealing 
with the single polymer diffusion) it is known since the work by Kirkwood and Riseman |23j . As to the internal poly- 
mer dynamics, our approach corresponds to that by Dubois- Violette and de Gennes [T(i| who, however, have assumed 
that the internal modes of the polymers should behave as the Zimm modes (i.e. with the dispersion ~ P~ 3 ^ 2 ), thus 
neglecting the Rouse contribution, initially being present in their theory. Such a simplification requires quantitative 
arguments and in many cases it is not substantiated, as well as the assumption of the continuous distribution of the 
internal modes; we believe that it is clearly shown in the present work. To take into account the presence of other 
coils in the solution, we have used the well-known Debye and Bueche (or Brinkman's) theory for a porous medium. 
Again, this approach has been already used in the polymer physics. However, coming from our earlier results on the 
hydrodynamic theory of the polymer dynamics, we could in summary build a model that is able to predict new results 
on the fundamental characteristics of polymer behavior in dilute solutions. Some of the quantities characterizing 
the polymer solutions (viscosity, Huggins coefficient) could be verified in standard experiments (with the necessary 
account for the draining parameter) . We have also proposed a description of the time dependence of the tendency to 
hydrodynamic screening in dilute polymer solutions; this effect seems to be suitable for computer simulations studies 
similar to those in Ref. Q- 
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